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On the Spectrum of Energy in Turbulent Shear Flow 

C. M. Tchen 

A turbulent flow with a given spatial pattern of mean motion invariant with time is 
treated. No solid boundary is considered. By means of the Fourier analysis of the Navier- 
Stokes equation, it is possible to obtain an equation of energy balance in spectral terms. 
This involves the following functions: Viscous dissipation, transfer, production, convection, 
and diffusion. The first two functions are common to the isotroj)ic and homogeneous turbu- 
lence without mean motion. Assuming that eddies controlled by inhomogeneity are much 
larger than those controlled by anisotropy, stationary spectral laws are obtained for small 
eddies in the viscous and nonviscous ranges. The agreement between theory and experiment 
is shown, 

1. Introductory Remarks and Hypotheses on the Energy Distribution 

In the statistical theories of isotropic turbulence, the irregular motion of the fluid is 
described by some characteristic functions: the correlation function, and the spectral function. 
The statistical theories devoted to the analysis of the correlation function were developed by 
von Karman, Burgers, and others [1].^ The study of the spectral function is based on the analysis 
of an irregular motion into harmonic components, as was developed by Burgers, for a model of 
turbulence [2]. The nonlinear process of the balance of turbulent energy is controlled by a 
transfer function, describing the transfer of energy from large eddies into small ones. Certain 
shapes of the transfer function were studied by Onsager [3], Obukhoff [4], von Weizacker [5], 
and Heisenberg [6], all leading to a spectral law in agreement with the Kolmogoroff [7] law of 
locally isotropic turbulence for the nonviscous range. It seems that the spectral law in isotropic 
turbulence is well determined, at least for the above range of small eddies. In the following 
pages an attempt is made to study the spectral law of small eddies in turbulent shear flow with 
a given pattern of mean motion without sohd boundaries. In addition to the viscous dissipation 
function and the transfer function, as occurring in isotropic turbulence without mean motion, 
we have to consider the production function, the convection function, and the diffusion function, 
which arise from the presence of a pattern of mean motion. 

The Navier-Stokes equation applied to a turbulent flow can be divided into two equations: 
an equation governing the mean, or primary motion C/^, and an equation governing the life 
history of the turbulent velocity, or velocity fluctuation u^ The equation for the turbulent 
velocity that interests us the most is, for an incompressible flow: 

-^+Wj^^'U,)=.-^-^^{Ua^)-^^{u,u,-u,u,+^^5,,). (1) 

The mean velocity Uf and the velocity fluctuation Uf must fulfill the equations of continuity 

bUilbXi=0, c)Ui/c>Xi = 0. (2) 

The independent variables in the equations are the time t and the coordinates Xi, Xj. Like the 
velocity, the pressure is decomposed into a mean pressure P, and a superposed pressure fluctua- 
tion p, p represents the density of the fluid and /x its viscosity, while v=ixlp. p, m, v are treated 
as constant. The indices i, j*, denote projection components. They are i=\, 2, 3,j = l. 2, 3. 
A summation is understood where the indices repeat. The bar over a quantity represents a 
spatial mean value. 

We shall study the spectral distribution of the energy of the turbulent motion in a given 
field of primary motion. The interaction between the primary and the turbulent motion is 
controlled by terms such as d{UiUj)ldxj, and d(UiUj)dXj in the hydrodynamical equation (1). 
According to this equation, the structure of the field at a given instant is characterized by three 
kinds of motions — a mean or primary motion, which gives rise to production, convection, and 
diffusion of turbulence; a turbulent motion, which possesses a transfer of energy ; and a molecular 
motion playing the role of viscous dissipation. The three motions have their own geometrical 
patterns generally different in average scales. In order to fix our ideas, we can choose a large 

1 Figures in brackets indicate the literature references at the end of this paper. 
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region of size R, in which the primary motion shows ceriain geometrical patterns. At a given 
point we draw an element of volume of dimensions 2Xi, 2X2, 2X3, having several times the 
dimensions of the largest eddy to be included in the study. If we study the spectrum, or the 
energy carried by the eddies of different sizes contained in the element of volume, we can 
expect, according to the differential orders of the terms in the hydrodynamical equation, that 
the molecular action enters 'from the end of the small eddies, and the action of the primary 
motion from the end of the large eddies. 

If we do not go down to the range of the molecular sizes, the action of the viscosity on 
the turbulent motion has a simple form, where the two motions are separated. On the other 
hand, if we do not go up to the very big eddies, and if the primary motion has a slowly varying 
pattern, or simple pattern, so that within the element of volume the velocity change can be 
considered constant and small, the action of the primary motion on the turbulent motion has 
also a simple form, where the two motions are separated without resonance. On the contrary, 
if the primary motion has a complexl}^ varying pattern, or rich pattern, so that in the volume 
element are presented various scales of the patterns of the primar}^ motion, and that certain 
small scales are comparable to the eddies of the turbulent motion, the two motions will interact 
with resonance and are not separable. Such circumstances will occur when the average scale 
of the primary motion, as depending on U^/U^\ U'^ jW" , etc., according to the similarity 
theory, is small and becomes comparable to the scales of the eddies to be studied. 

Let us isolate the finite element of volume from its surroundings and study first, by means 
of an elementary method, the balance of energy per unit volume among the different eddies. 
For this purpose, let k be the wave-number magnitude, v^ the velocity, and l^l/k the mixing 
length. The vorticity of the turbulent motion can be represented by vJL The vorticity 
presented in the pattern of the primar}- motion, averaged over the element, is supposed to be 
given, and is denoted by U\ 

Let Sjc be the flow of energy from larger eddies (larger than the eddy k) to smaller eddies 
(smaller than the eddy k). Such a flow of energy corresponds to a dissipation due to the 
following exchanges of momentum: between the molecular motion and the turbulent motion, 
between the small and large eddies, and between the turbulent motion and the primary motion. 
According to Prandtl and Boussinesq, it is known that the dissipation is the product of the viscos- 
ity with the square of the vorticity. 

For the dissipation of molecular origin, we can write ijlvI/P, where fi is the coefficient of 
viscosity due to molecular agitations. For the dissipation of turbulent origin, that is, the 
dissipation due to the exchange of momentum between the eddies of the turbulent motion 
itself, and between the turbulent motion and the primary motion, we can write 

turbulent dissipation^p^^:- compound square velocity. 

Here v^ is the turbulent viscosity of Heisenberg Vk=Vjcl. If there is no resonance between the 
primary motion and the turbulent motion, the compound square vorticity is formed by the 
following independent terms: 

the interaction term (vjl) U' being dropped out. If there is strong resonance between the two 
motions, the independent term C/^^ must be dropped out, because any primary vorticity field 
that is able to create momentum exchange with eddies must be coupled to the turbulent 
vorticity. Thus the compound square vorticity becomes in this case 



Hence the equation of total dissipation is 



Sic=li -p+P^A;0|| + < 



/U' \ no resonance. 
^-^ U'\ strong resonance. 



The case of intermediate resonance is not investigated here. 

A spectral function F{k) can be introduced such that pF{k)dk is the turbulent energy between 
k and k-\-dk. Then by means of dimensional reasoning, the dissipation equation can be put 
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ill spectral terms as follows: 



^=2^ (^dk/k''F(k') + 2K rd¥\K^ rdk'k'^F{¥) + 
P Jo Jk V k'^ Jo 

\-7 2 a A: A: r(A: ) strong resonance, 

where k, k\ k" are numerical constants to be determined by experiments. 

We shall study the spectral laws in the equilibrium range of spectrum, which, as defined 
by Heisenberg, is the range controlled by a constant flow of energy Sk, independent of k. In the 
equilibrium range, two subranges are defined: the nonviscous, or inertial, subrange, as charac- 
terized by Vjc/kv :^l, and the viscous subrange, as characterized by Vk/kv <^1. For the non- 
viscous subrange we obtain 



and 



V ^k-"% 


or 


F^k-'i^ 


Z);-.'^ constant 




or F'-^k'^ 


ous subran<'e 


we 


obtahi 



for small U^ and no resonance; 

for lar^-e r^ and stronti' resonance. 



ih'^k ^, or F^k ^ for arl)itrary f/' and resonance. 

The spectral law k~^ is not expected to go down to molecular scales, as th(»n the strong inter- 
ference between the mohn-ular motion and the turbulent motion will nhcv the al)ove simple 
forms of the molecidar and turbulent dissipations. 

The above elementary study of the spectrum of turbulence is based on the considerations 
3f the different forms of the exchange of momentum within an isolated volume element. Actu- 
ally, the volume element is not isolated, and the general character of the fi(4d from one element 
to another may vary. Under- those circumstances we may expect that the gc^uM'al sp(H'trum 
depends on the position of tlie eleuK^nt analyzed, and that there exists an exchange^ of momen- 
tum between the (dement chosen and its neighbors, in thc^ forms of conviH'tion and diffusion. 
However, the kit ter forms of excliange, called inh()m()^(Mi(H)us (\\chang{\ hav(^ scales c()mj)arable 
to the dimensions of the volume elements and thus affect predominantly the large echlies 
rather than the eddies in the equilibrium range in which we are interested for the present. 
We shall return to the discussions of convection and diffusion in sections 3 and 4. We con- 
clude that in the study of the equilibrium range of spectrum, we can analyze the inhomo- 
geneous field as follows: 

1. We choose a finite element of volume 2Xi-2.Y2-2A'3, which is of the order of several 
times the scales of the lai'gest edd}^ of the equilibrium range. In the element of volume a 
turbulent motion is suj)erposed to a primary motion exhibiting certain geometrical ])at terns. 
Both motions can be analyzed into Fourier components by the same method as a{)pli(Hl to 
a homogeneous field. 

2. The inliomogeneity of the general field creates an inhomogeneous exchange of momen- 
tum between the chosen element of volume and its neighbors, resulting in convection and 
diffusion. Such an exchange of large scales is not important to the study of eddies in the 
equilibrium range (small eddies). Thus it can be estimat(^d approximately and added to the 
homogeneous exchange as small corrections. 

2. Harmonic Analysis of a Turbulent Motion 

The results-predicted in section 1 by means of elementary considerations of energy distribu- 
tion can be obtained by a harmonic analysis of the hydrodynamical equations. Such an 
analysis will reveal some important details of the exchange mechanism and of the energy 
balance, which are rather obscure in the elementary considerations. 

The applications of Fourier analysis to a turbulent field has the difficulty that Fourier 
integrals can be defined only for f mictions that vanish at infinit}^ in such a way that their squares 
take a finite integrable value. To overcome this difficulty, one usually assumes in the harmonic 
analysis of turbulence that we can restrict to the values of the function within a finite domain, 
say a finite element of volume 2A"i-2X2-2X3 taken at a point Xi, X2, Xz as center, while outside 
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this domain the function is replaced by zero.^ For the purpose of the analysis of turbulence, 
the Fourier integral of the turbulent velocity in such a finite domain will be sufficiently represen- 
tative of the energy distribution of eddies smaller than the element. The amplitude functions 
of the Fourier coefficients are then obtained from 



I rx+x , 

{ZTyjx-x 



(3) 



1 r^+^ 



-ik-x 



Ai{k)=^A c?x'C/i(Jc')e-'"-^ (3a) 

c(k) = j^, r^^'c/jcXjcOe-'"-^' (3b) 

{^-n-) Jx-x 

with the equations of continuity 

kMk)=0, k,A,{k)^0. (4) 

Here Ui, Ui, p are real quantities representing the turbulent velocity components, the mean 
velocity components, and the turbulent pressure, respectively, af(fe), Ai(k), c{k) are the corre- 
sponding Fourier coefficients at the wave-number vector k ( A:= 2 tt /wavelength), and are complex 
quantities, with a*{k)=ai(—k)y etc. An asterisk over a quantity denotes its conjugate, and 
t=^J—l. The following symbol for the volume integral is used: 



I dx^= dx[dx2dxs. 



The actual change of Uf and p in the course of time is resolved into the change of the Fourier 
coefficients ai and c with the time. Provisionally, we restrict to the singie instant. The 
analysis in a finite volume renders the Fourier coefficients bounded. Outside and on the 
boundary of the finite volume {x^<x—X, jk:'>jc+^), the function Ui{x^) and its derivative 
bui{x^)l'dx] are supposed to vanish; so also do p{x') and Ui{x'), Such functions are called 
truncated junctions. 



To show the significance of the limit X, we form: 



Ui{x)Uj{x+r)= r^dk' {'^^dk''ai{k')aj{-k''\ 

J — 00 J —CO 



y 



,i{k'—k")'X—i,k"'r 



The mean value Ui(x)Uj{x-\-r) is obtained by integrating with respect to JC from —X' to -\-X^ , 
where X'>X. We have 

Ui{x)u^{x-\-r)= ^ i ^ r^dk' r^ dVQ-^^"''ai{k')a^{-k')' 
(sin k[ — k"^X[ sin (fcg— ^2)^2 sin {y^—¥^X'^ 

h' — V h' —V h' — h" ' 

which for Jf'— >oo leads to 

Ui{x)Uj{x+r)=j^YY dkai{k)aji—k)e-'^'\ 



Since Ui{x)uj(x+r) must be independent of X, it follows that the absolute value of the ampli- 
tude function must be proportional to {XiXoX^y^^. As at is obtained by the summation of 
a large number of variables like Ui{x')e'^'^' according to (3), the fact that its absolute value is 
proportional to the square root of the volume must be related to the same results found in 
problems of radom walk, even though the paths may be correlated [10]. 

2 To avoid discontinuities that may arise by such a method at the edges of the elementary volume, we may instead multiply the function by 
a certain function that is unity within the volume, and decays rapidly to zero beyond the volume. 
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If we take r=0 and j=i, we find 

^i = X XoX I- ^^ cii{k)ai{—k). 

Instead of extending the integral over the whole fc-space, we can integrate over a spherical 
shell of radii between k and k+dk, and we find 

2F{k)= ^ !L' r^dcp Cdd sin ek'ai{k)aii-k) with 2 rdk'F{k')=u^ (5) 

A1A2A3J0 Jo Jo 

Here k^sinOdddip is the volume element with polar angles 6, (p, and F(k)dk is the energy associated 
with the wave numbers ranged between k and k-\-dk. F{k) will be referred to customarily 
as spectral function of turbulent energy. Also we can define a spectral function of the turbulent 
shear stress Fij{k) such that 

2Fij(k)= ^ "C X. C'^dip rdesmekmai(k)aj(-k) + ai(-k)aj(k)] with 2 (^ dk'Fij{k')=v~Uj. 

A1A2A3J0 Jo Jo 

(6) 
By a Fourier transform of (1) term by term, we have 



— ^7-^+ I dniniai{n)Aj{k—n) = — vk^ai{k)— j dnLnjAi(n)aj{k—n) — 

J 4-00 ^ 
dninjai{n)aj{k — n) ttiC(A:) — inhomogeneous terms, 
p 

where the inhomogeneous terms arise from the inhomogeneous exchange as explained in section 
1 and to be determined in sections 3 and 4. In the above equation c can be eliminated by its 
value obtained from the above equation multiplied by ki, subjected to the conditions of con- 
tinuity (4). The transformed equation is then 

^^^=-vk'a,{k)- r"c/nt%a/fe-n)ra,(n)-^a,(n)1- 

j dnLnjaj{k — n)\ Ai{n) — ^As(n) — 

J +00 r- k k '1 . 

dnLnjAj{k — n)\ ai{n) — j-^as{n) —inhomogeneous terms. 

An equation describing the rate of change of the energy distribution can be obtained by forming 
the differential quotient {d/dt)ai{k)ai{—k) as follows: 

■^ ai{k)ai{— k) = — 2v¥ai{k)ai{— k) — 

dnirijaiin) [aj{k—n)ai{—k) + aj{—k—n)ai{k)] — 

J — 00 

dnLnjAi{n)[aj{k—n)ai(—k) + aj{—k—n)ai{k)]— (7) 

J —CO 

dmnjaiin) [Aj{k — n)ai{—k)+Aj{—k — n)ai{k)]—\v^ioin.oge>nQOus terms. 
»/ — 00 

The physical meaning of the equation for energy distribution (7) can be compared with 
the following equation for over-all balance obtained by multiplying (1) by Ui and taking the 
spatial mean values as follows: 



dt 2 
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Equation (7) represents the balance of energy associated witli the wave number k, while eq (7a) 
represents the balance of energy covering the whole spectrum. The correspondence between 
terms of the two equations is as follows: 
(a) Rate of change of energy: 



I dk-^^ai{k)ai{—k) 



2 X,X2XJ-^ ^^-^V'V-nV 'V ^^ 2 



For a stationary flow these two terms vanish. 
(b) Dissipation by viscosity: 



2 ^iA2^3j— CO 

(c) Transfer of energy from eddies to eddies: 



^ f ^ dk2vk^ai{k)ai{-k) = v{bUilbx,y. 



-J 



■+CC 

dn L njaiin) [a^fc — n)ai{— k) + aj{—k — n)a^(fc)] . 



This term which figures in (7) has no corresponding term in (7a), because the total transfer 
must vanish for the whole spectrum. 

(d) Production of turbulent energy by shearing stresses: 

t; V V V dJk dnLnjAiin)[aj{k — n)ai{—k) + aj{—k—n)ai{k)] = UiUj-^^' 

(e) Convection of turbulent energy by interaction with mean motion: 

2 X X X J- ^^ _ ^^^^''^i^^(^)[^i<^-^)^^(-*^)+^X-fc-^)«z(^)]=^^.(^ ^y 

(f) Inhomogeneous transport: This is the part of convection from the inhomogeneity 
of the field. Thus in (7a) the terms 

correspond to the inhomogeneous terms of (7), to be determined in sections 3 and 4. For a 
stationary condition the left-hand side of (7) must vanish. Integration between k and oo 
gives 

2p {'^dk'k'^F{k') = 2v [ \lk'¥~F{¥) + W{k, oo) + xi,(k,oo)+C{k, oo)+K(A:, oo). (8) 

Jo Jo 

This equation expresses the ecpiilibrium in the exchange of energy. The left-hand side is 
the total viscous dissipation covering the whole spectrum. It is the total flow of energy. 
It is independent of k and will be represented by e. The first term on the right-hand side 
gives the loss of energy through viscous dissipation associated with wave numbers not ex- 
ceeding k. The other four terms give the flow of energy from this part of the spectrum through 
interaction with the rem.ainder of the spectrum., respectfully, in the formes of transfer, produc- 
tion, convection, and diffusion, as to be obtained by integrating the corresponding terms of 
(7) for wave num.bers not exceeding k. 

3. Transfer, Production, and Convection of Turbulence 

For the transfer function W{k, oo)^ we assume the Heisenberg formula [6]. 



W{k, ^) = Kf^dk"^^ 2£dk'k''Fik'), (9) 
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where /c is a numerical coiisiaiit. The justification of the dissipative cliaracter of the transfer 
and hence of the Heisenberg theory, can be sought from the transport processes by turbulence 
on the basis of the hydrodynamical equations. We shall return to this problem in a subsequent 
comroiniication. 

For the production function ^(/:, oo), we assume a gradient of Ux in the 2*2 direction, as 
expressed by bUilc^x-z, all other gradients being negligible, and distinguish the case of no res- 
onance and strong resonance, depending on the scale of the prim.ary motion, as mentioned in 
section 1. 

(a) When there is negligible resonance between the mean and tbe turbulent motions, we 
can neglect n compared to k in the following production term: 



X 



+ 00 

dnLn2^i(n)[a2{k — n)ax{—k) + a2(—k'-n)ai{k)]^[ai{~k)a2{k) + ai{k)a2{--k)] -U^ 



where ^ 

^+co 

dnL7i2Ai(n)=U' 



f- 



is the magnitude of the average vorticity of the mean motion for a quickly (h*ot)])in.g spectrum 
Ai{n). Integration of the production term over the wave numbers of magnitudes between 
k and 00 gives 

^{k,oo) = 2 rdk'Fr2kk')'U\ (10) 

If, according to l^randtl and Broiissiiiesfj [cS, 9,], we wi'ite 

/»00 , 

where k' dk"SF{}c")\k"^ is a turbulent viscositv, and k is a numerical constant, ecj (10) is re- 

r°° \fW^ 2 

^(Z:-, oo):=r /<' ^/^''^_A__^ . u' ^ for negligible resonance. (12) 



duced to 



(b) When there is strong resonance between the mean and turbulent motions as described 
in section 1, in such a way that we cannot separate them a priori by the method used in (a), 
a new m.ethod has to be devised in order to formulate the production function ^. Integra- 
tion of (7) outside the sphere of radius k gives the transfer and ])roduction functions as follows: 

^3 /"ex, /»co 

— W{k,co) = -—--—^ dk'' (/ntn/tj(n)a,(A:'' — n)a,( — fe'0 + <^onjugate 

AiAoAsJit JO 

3 /^oo /»00 

— "^{k, 00)=: dk'' dn injA,{n)aj{k'' — n)a^(-~k'') + conjugs.te. 

Ji]A2A^J k Jo 



2 By applying the Fourier formula (3a), the amplitude function of the vorticity is 

1712^1(712) = — dX2- ^ ^ 



the factor before (7' being unity. 



1 r+''^^ 

•''J-X2 
J_^ rf^^2C 7*2-1.(7^2)=- J_^ ^''U-X2 ^'^^~^' 

vithin (—X2,-\-X2), and has an average value U', as if 
I (y7i2t 7?2-4i(r?2)= — I dn2 U'=U\ 

J-co ^J-co n2 



and. 



If the vorticity does not change too mucli within (--X2,-\-X2), and has an average value U', as is the case here, we can write, as a first approxi- 
mation: 
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We notice that they are analogous in structural form, except that the transfer function contains 
a triple correlation 

q{n, fc") — ^% ai(n)a^(A:"— n)ai(— fe'0 + conjugate, 

which is antisymmetrical: qi{n,k") = —g{k",n), while the production function contains a 
triple correlation involving the gradient of the mean motion and is not antisymmetrical. For 
the antisymmetrical function, we can replace 

dk" dn by dk'' dn, (13) 

Jk Jo Jk Jo 

in agreement with the form of integrals given to W(k, co) in (9). For the production func- 
tion, we can take the following form, because of its structual analogy with the transfer function: 

dk^^-xj ^^3 • (vorticity of turbulent motion • vorticity of mean motion) (14) 

ft ^ tC 

Now the vorticity of the turbulent motion can be written in the usual form 

[2rdk'k''F{k')'\% 

while the vorticity of the mean motion is U\ We can remark that transformation (13) does 
not apply to correlations involving the mean motion, and therefore the vorticity of the mean 
motion is written for its whole spectrum. After substitution for the vorticities, (14) becomes: 

<if(k, oo) = k" r dk\l'^^[2 r dk'k''FikyU''f' (15) 

Jk V /:" L J/c J 

Finally we shall integrate the convection term of (7) for the region outside the sphere of radius 
k. We can remark that the transfer, production, and convection contain, respectively, the 
triple correlations 

Lnjai{n)aj(k — n)ai{—k)-\-conjvigSi,te 

LnjAi{n)aj(k — n)ai( — k) + con]ug8iie 

Lnjai{n)Aj(k — n)ai(—k) + conjugate. 

By analogy with the expression (9) and (15) for the transfer and the production, we can for- 
mulate an expression for the convection function. However, we assume that the apparent 
viscosity is much sm.aller for the m.ean motion than for the turbulent motion, for the sufficiently 
large value of k herewith considered; hence the convection function is negligible compared 
to W{k, 00 ) for sufficiently large k. 

4. Diffusion by Inhomogeneity of the Field 

Let us consider a case with the following simple conditions: 

(a) The diffusion due to the inhomogeneity is primarily a large scale phenomenon, and 
hence the molecular diffusion can be neglected compared to the turbulent diffusion. 

(b) The mean pressure is constant in time and space. 

The difference of energy distribution in the space causes a flux of energy as a consequence 
of the exchange of matter. This flux of energy can be written as 

flux of enevgy=r' j^dk"^^ • ^ [p pZ:'/'(fc')]- 

The first integral is the turbulent viscosity, the gradient is the gradient of the energy distribu- 
tion, and /c"' is a numerical constant to be determined by experiment. By taking the differ- 
ence of the fluxes across opposite faces of the elementary space and dividing by p, we obtain the 
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rate of decrease of energy by diffusion per unit mass and unit time equal to 



Finally, if k is chosen sufficiently large so that the range from to k includes practically all 
the energy gradient, we can write 

^rdk'F{k')^^(^^\ 
OX2J0 0X2\2 V 

Hence 



«'.")=". 41- j:-Vf? ^(1=3)] 



5. Energy Spectrum of Turbulence 



(16) 



In the equilibrium range the energy spectrum is determined by eo 8, which becomes, 
after substitution of (9), (12), (15), and (16), 

e=2v r d¥¥'F{k')+2K r dk'\i^^ r dk'k''F{¥)+ 

Jo Jk V k" Jo 



with 



f=< 



1, no resonance 




k" 




i 


7 


) strong resonance 



(17) 



(17a) 



We shall study the energy spectrum for the following subranges: (a) non viscous subrange 
with large U^ and strong resonance, (b) nonviscous subrange Avitli small U^ and no resonance, 
and (c) viscous subrange with arbitrary U^ and arbitrary resonance. 

(a) In the nonviscous subrange of the spectrum, if there is a strong resonance and if U' 
is large, along with sufficiently small i, we have 



I Vdk'¥'F{¥)<:U'\ 
Jo 



thus the production function is predominant, and all other functions are negligible. In this 
way eq (17) reduces to 



AT'=i?'"^Wb^l''''"^^''^l 



and has the solution F={elK"U')k~^ to be obtained simply by introducing the new variable 
fr=2\ dk'k'^Fik'), If the diffusion term is small, but not negligible, then as a first approxi- 
mation, the solution is 



F = ^^,k-\l-krlk), 



where 



h- 



k" d 



e ^x\\^"UV ^2-2(2 ""vj 



and krjk is assumed to be small compared to unity. 
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(b) In the non viscous subrange, if there is no resonance and U' is small, along with suffi- 
cientl}" large k, we have 



2 dk'¥'F{k'):^U'\ 



Thus the production function and the diffusion function are small, although not negligible 
com.pared to the transfer function, and eq (17) gives the following solution as a first approxi- 
mation: 

F=(^^"k -^/^ [1 - {k,ikyi'\ ''', (1 8) 

vN^nere 

,,.4^-.)V+(|^i)'"A[..,.^^(i^)]. 

ana k-pjk is assumed to be small compared to unity. The spectrum 

ox tlie homogeneous and isotropic turbulence is a particular case corresponding to ^^ = 0. 

For the reasons incorporated already in the conditions under which the respective laws 
k~^ and k~^'^ are obtained, it is to be understood that the former law refers usually to larger 
eddies than does the latter one. Depending upon specific conditions, a spectrum may show 
clearly a portion representative of one of the laws and not the other, or ma}^ show clearl}- two 
portions representative of both laws successively. In the latter case, the primary motion 
possesses a pattern in resonance with the secondary motion. As the result of the resonance, 
the production function is first larger than the transfer function at small k, hence the spectrum 
shows a k~^ law. Afterwards, as k increases, the production function becomes smaller than 
the transfer function, while U' is not large enough to com.pensate, hence the spectrum changes 
over to a k~^'^ law. The ranges of validity of the laws depend upon the ranges of predominance 
either by the production function or by the transfer function. For the viscous subrange, which 
we shall study in (c), the eddies are m.uch sm.aller than the nonviscous subrange. Consider the 
case where we have a i"^ law, thus k is sm.all and U' is large. Suppose that the Reynolds 
num.ber is not sufficiently large, so that the viscous dissipation m.ay becom.e so im.portant as 
to alter the tail portion of the k~^ spectrum. Under such circumstances, for increasing wave 
nmnbers, the k '^ law may pass over directly to the k"'^ law, without a noticeable portion of 
the k^^'^ law. For very sm.all Reynolds num.bers, a larger part of the tail portion or even the 
whole of the k~^ law m.ay be affected, as could be observed in the turbulent flow in the close 
vicinity of a wall. 

(c) In the viscous subrange, the turbulent friction of turbulent viscosity 



K {^dk\ 

J k 



is much smaller than the molecular friction. By differentiating (17) with respect to k, and 
b}^ neglecting the turbulent viscosity com.pared to the m.olecular viscosity, we have 

V \ k^l_ K e K e dxl^ ' _ 

where f^ is obtained by putting k= oo in expression (17a). 
The solution follows: 

For homogeneous and isotropic turbulence, the above formula degenerates to the Heisenberg 
formula [6] 

60 



v^ 


•t> 










\ 












\»\ -1 

\ * 










o 


• 
• 










o 

• 










\ 4 















Figure 1. Energy spectr^um in a turbulent boundary layer. 

Fn(ki) is the spactrum of uy^ associated with the wave number ki. The free stream velocity is 15 m/sec. The measurements are made at a distance of 320 cm 
downstream from the leading edge of the flat plate. The thickness of the boundary layer is 7.6 cm. The open circles are measurements at a distance of 6 cm. from 
the wall, and the cbsed circles are measurements at a distance of 0.4 cm from the wall. The data are new measurements made by KlobanotT at the National Bureau 
of Standards. Some early results are published in reference [11]. 

The /j-"^ law follows from the assumption of the Heisenberg expression for the transfer function. 
We see that the eddies in the viscous subrange of the spectrum, are not effected by any shear. 

We shall compare (see fig. 1) the above results with the experim.ental investigations in a 
turbulent boundary layer m.ade by P. S. Klebanoff and Z. W. Dielil_^[ll]. It is to be rem.arked 
that the available measurements give the spectral function for v^/2 in the one-dim.ensional 
analysis of the wave-num.ber space, while the spectral function in the theory refers to the 
energy u?/2 for the wave-num/ber magnitude in the three dim.ensional analysis of the wave- 
number space. Assuming that the power laws for sufficiently large k rem.ains unaltered by the 
transform.ation between the two kinds of spectral functions for the cases considered, we see 
that the laws k~^, k~^^^, k~'^ agree with experim.ents respectively for the following cases: 

(a) Near the wall where the average scale of the primary motion is small, there is strong 
resonance between the primary motion and the turbulent motion. If, moreover, 



t7''>2 Cd¥¥^F{k') 



in the non viscous subrange, the production of turbulence predominates, and the spectral law 
is k~^. 

(b) Far away from the wall, and still not close to the rough edge of the boundary layer, the 
average scale of the primary motion is large, and the resonance is negligible. If, moreover, 



t/''«:2 Cdk'k''F{¥\ 



in the nonviscous subrange, the transfer of energy among the eddies predominates, and the spec- 
tral law is k~^'^. 
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(c) In tlie viscous subrange, the eddies are under no circumstance effected by the shear. 
Therefore, everywhere in the boundary layer the power law k~'^ of Heisenberg is reproduced, 
although this law is not expected to continue to infinite wave numbers. 

If we consider as before a turbulent flow with a mean velocity in the Xi direction, with a 
gradient d?7i/d^2, and write energy equations of the type (7) and (7a), we notice that the 
production function will be absent in the energy equations for the U2 and Uz components. 
Therefore, the spectra of U2 and Uz will not show any k~^ law which is characteristic of the 
production function. 
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